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We have carried out a neutron scattering study of the instantaneous spin-spin correlations in 
La2CuC>4 (Tjv = 325 K) over the temperature range 337 K to 824 K. Incident neutron energies 
varying from 14.7 meV to 115 meV have been employed in order to guarantee that the energy inte- 
gration is carried out properly. The results so-obtained for the spin correlation length as a function 
of temperature when expressed in reduced units agree quantitatively both with previous results for 
the two dimensional (2D) tetragonal material Sr2Cu02Cl2 and with quantum Monte Carlo results 
for the nearest neighbor square lattice S=l/2 Heisenberg model. All of the experimental and numer- 
ical results for the correlation length are well described without any adjustable parameters by the 
behavior predicted for the quantum non-linear sigma model in the low temperature renormalized 
classical regime. The amplitude, on the other hand, deviates subtly from the predicted low tem- 
perature behavior. These results are discussed in the context of recent theory for the 2D quantum 
Heisenberg model. 
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The physics of low dimensional quantum Heisenberg 
antiferromagnets has been the subject of research ever 
since the advent of modern quantum and statistical me- 
chanics [1,2]. Interest in two dimensional (2D) systems 
was heightened by the discovery of high temperature su- 
perconductivity in the lamellar copper oxides [3] . Specif- 
ically, it was realized early on that the parent compounds 
such as La2Cu04 correspond to rather good approxima- 
tions to the S=l/2 2D square-lattice quantum Heisen- 
berg antiferromagnet (2DSLQHA)[4,5]. It seems at least 
possible that the 2D magnetism may in some way be 
essential to the superconductivity in the charge carrier 
doped cuprates. Further, the magnetism itself is of fun- 
damental interest as a quantum many body phenomenon 
in lower dimensions. 

Early experiments by Endoh and co-workers [5] showed 
that over a wide range of temperatures above the three 
dimensional Neel ordering transition in L^CuC^+j, (that 
is, La2Cu04 with a small amount of excess oxygen) the 
instantaneous spin-spin correlations were purely two di- 
mensional and that the correlation length diverged expo- 
nentially in 1/T. This led to a flurry of theoretical activ- 
ity [2] including most especially work based on the quan- 
tum non-linear sigma model (QNLerM) by Chakravarty, 
Halperin and Nelson (CHN) [6] and Hasenfratz and Nie- 
dermayer (HN) [7]. These theories are all based on the 
2D Heisenberg Hamiltonian which for nearest neighbor 
(nn) interactions alone takes the form 

H = J Si'Si+Kn (1) 

<i,5 nn > 



where the summation is over nn pairs on a square lattice. 

In La2Cu04, for temperatures below the tetragonal 
(lA/mmm) - orthorhombic (Bmab) structural phase tran- 
sition temperature of T st = 530 K, the leading terms in 
the spin Hamiltonian [8,9] are 

H = J ( Si ■ Si+§ nn + OLnnn /J Si ■ Si+$ nnn 

<i,5„„> <i,$nnn> 

+ a xy Y S i S i+S nn + a± J Si ■ S l+5±j 

<i,S nn > <i,8 ±j > 

+ ttDM (-) l a-S z x S i+5nn ). (2) 

Here, a nnn , a xy , a±j, and aoM represent the re- 
duced next nearest neighbor in-plane Heisenberg ex- 
change coupling, XY anisotropy, interlayer coupling and 
Dzyaloshinski-Moriya antisymmetric exchange, respec- 
tively, and Sf is the c component of the spin at site i. 
The fourth term in Eq. (2) explicitly includes the two dif- 
ferent out-of-plane neighbors at 5±% and <5_L2. Note that, 
as was implicit in the work of Thio et al. [8] , the sign of 
the antisymmetric term changes on opposite sublattices 
because of the opposite rotation of the CuOg octahedra. 
This Dzyaloshinski-Moriya term originates from a small 
rotation of the CuOe octahedra about the a axis. In 
the tetragonal phase cxdm = and the nearest neighbor 
out-of-plane effective coupling vanishes since an = aj_2. 

The most complete experimental study to-date is on 
the material Sr2Cu02Cl2 [10] rather than La2Cu04. The 
reasons for this are twofold: First, Sr 2 Cu02Cl2 is very 
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difficult to dope so that there are no complications aris- 
ing from the effects of doped electrons or holes on the 
spin correlations. Second, since S^CuC^C^ is tetrago- 
nal down to the lowest temperatures measured (< 10 K), 
c*dm=0 and the nearest neighbor intcrplanar coupling 
vanishes to leading order, that is a±\ = a± 2 - As shown 
in Table (1) there is a small XY anisotropy. In addition, 
from results in Sr 2 Cu304Cl2 [11], we infer that there is a 
next nearest neighbor in-plane Heiscnberg exchange cou- 
pling which is about 8% of the nearest neighbor value. 
To first order, the latter should simply lead to a slight 
renormalization of the effective J in Eq. (1). The XY 
anisotropy will lead to a crossover from Heiscnberg to 
XY behavior for correlation lengths £/a ~ 100. Thus, 
Sr 2 Cu02Ci2 should be a good realization of the S=l/2 
2DSLQHA for length scales ~ 100. This has, in fact, been 
confirmed in detail experimentally [10]; specifically, over 
a wide range of length scales the 2D correlation length 
measured in Sr2Cu02Cl2 agrees quantitatively with re- 
sults from quantum Monte Carlo (QMC) calculations 
carried out on the Hamiltonian Eq. (1) with S=l/2 [12- 
14]. The value for J for Sr 2 Cu0 2 Cl 2 listed in Table (1) 
is deduced from two magnon Raman scattering measure- 
ments [15]. 

Both the QMC and the Sr 2 Cu0 2 Cl 2 experimental re- 
sults for the correlation length in turn are quantitatively 
predicted by theory based on the QNLaM in the low 
temperature renormalized classical (RC) regime [6, 7]. 
This comparison again involves no adjustable parame- 
ters. Surprisingly, this agreement holds for correlation 
lengths as short as a few lattice constants. This is far 
outside of the temperature range where the QNLctM-RC 
theory should hold. A plausible explanation for this un- 
expected agreement has been given by Beard et al.[13] 

In spite of the fact that the progenitor of this work 
was the discovery of high temperature superconductiv- 
ity in La 2 _ 2: Ba 2: Cu04 [3] together with the early work 
on the 2D spin correlations in La 2 Cu04 +a [5], our 
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TABLE I. Neel temperature, superexchange energy, and 
corrections to the 2D Heisenberg Hamiltonian for La2CuC>4 
[9] and S^CuChCL; [10]. oldm and olxy are larger than the 
values quoted in references [9] and [10] by factors of (Z c /Z g ) 
and (Z c /Z g f respectively. Here Z c (l/2) ~ 1.17 and Z g (l/2) 
~ 0.6 are the quantum renormalization factors for the spin 
wave velocity and spin wave gap respectively. 



knowledge of the spin correlations in stoichiometric 
La 2 CuC"4 is rather limited. The primary correlation 
length data for La 2 CuC>4 originate from the neutron scat- 
tering study of Keimer et al. [9] on a carrier-free single 
crystal of La2Cu04 with Tn = 325 K. The Keimer et 
al. [9] data on the correlation length and structure factor 
extend up to 550 K. Their measurements are generally 
consistent with the Sr 2 Cu0 2 Cl 2 , QMC and QNLctM-RC 
results, but there appear to be systematic discrepancies 
at the limit of the error bars for the correlation length at 
both low and high temperatures. These neutron exper- 
iments were carried out using a single incident neutron 
energy of 31 meV. It seems likely that the discrepan- 
cies are an experimental artifact originating from the use 
of a single incident neutron energy over a wide range 
of temperatures. Alternatively, they could represent a 
real effect originating from the antisymmetric exchange 
and interplanar coupling terms in Eq. (2) for La 2 Cu04. 
Clearly, therefore, it is important to carry out a more 
complete study of the spin-spin correlations in La2Cu04 
in order to characterize fully the magnetism in this par- 
ent compound of the monolayer high temperature super- 
conductors. Such data would also be valuable for the 
interpretation of NQR results in La 2 Cu04 [16]. Finally, 
there have been some important advances in our under- 
standing of the theory for the 2DSLQHA since the work 
of Greven et al. [10] on Sr 2 Cu02Cl2 and it is therefore of 
value to re-examine the relationships between the results 
of experiments in real systems and theory. 

The experiments were carried out primarily on the H7 
triple-axis spectrometer at the High Flux Beam Reac- 
tor at Brookhaven National Laboratory. The measure- 
ments utilized the same single crystal of La2Cu04 as em- 
ployed by Keimer et al. [9]; this crystal had a volume of 
about 1.5 cm 3 . Throughout this paper we use Bmab or- 
thorhombic axes; at Tjv =325 K the lattice constants are 
a = 5.338 A, b = 5.406 A, and c = 13.141 A. We show in 
Fig. 1 the temperature dependence of the (0 1 2) nuclear 
supcrlattice peak intensity together with the reduced or- 
thorhombic splitting (b-a)/(b+a) [17]. As is evident from 
Fig. 1, the sample of La 2 Cu04 shows a sharp tetragonal- 
orthorhombic structural phase transition at T st — 530. 5± 
0.5 K. The sharpness of the transition in turn reflects the 
microscopic homogeneity of this sample. 

The magnetic neutron scattering experiments were car- 
ried out in the energy-integrating two axis mode. For 2D 
systems the integration over energy is carried out au- 
tomatically in a two axis experiment provided that the 
outgoing neutron wave vector kf is perpendicular to the 
2D planes and provided that the neutron energy is sig- 
nificantly larger than the characteristic energy oj of the 
spin fluctuations at a given temperature [18]. From the 
theory of CHN [6, 19] one has 
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FIG. 1. Orthorhombic splitting and (0,1,2) superlattice 
peak intensity versus temperature; the data are normalized 
relative to each other in the temperature region of overlap. 
The solid line is the result of a fit to a power law A(T st -T) 2 ' 3 
with (3 = 0.31 ± 0.01 and T st = 530.5±0.5 K. 

where c and p s are the zero temperature spin wave ve- 
locity and spin stiffness respectively. For La2CuC>4 this 
becomes [20] 
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Quantum Monte Carlo calculations by Makivic and 
Jarrell [21] at intermediate temperatures generally are 
well described by Eq. (4) but with an amplitude that 
is approximately twice as large. Specifically, Makivic 
and Jarrell [21] find that between 550 K and 800 K for 
J = 135 meV, as in La2Cu04 [20], ujq varies from <~ 
25 meV to ~ 63 meV. Accordingly, the following protocol 
was used in our measurements. Neutrons with incident 
energies Ej = 14.7 meV were used at the lower tempera- 
tures. With increasing temperature and hence decreasing 
£ the incoming neutron energy was progressively raised 
to 41 meV, 90 meV and 115 meV. To ensure that the 
energy integration was carried out correctly, it was re- 
quired that the results for the correlation length in the 
temperature regions of overlap agreed with each other to 
well within the experimental errors. 

We show first in Fig. 2 preparatory data taken at a 
temperature of T = 328 K which is just above the 3D 
Neel temperature of 325 K. The incident neutron energy 
was Ej = 3.6 meV which results in very high momentum 
resolution. The two peaks evident in Fig. 2 orginatc from 
the two rods of scattering which are along (1,0, I) and 
(0, 1,1), respectively. The equi-intensity of the two peaks 
implies that at 328 K the 2D spin fluctuations have at 
least XY symmetry, that is, at 328 K there is no measur- 
able in-plane anisotropy induced by the antisymmetric 
exchange terms in Eq. (2). 
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FIG. 2. E = 3.6 meV two-axis scan across the 2D rods 
at (1,0, 1.59) and (0, 1, 1.59). The solid line is the result of 
a fit to two Lorentzians, Eq. (5), centered about (1,0, 1) and 
(0, 1,2) respectively, convolved with the instrumental resolu- 
tion function. The fit gives £ _1 = 0.0011 ±0.0004 reciprocal 
lattice units. 

In Figs. 3 and 4 some representative energy-integrating 
scans for Ej = 41 meV and Ej = 115 meV are shown. 
The collimations were set to 20' - lO'-S-lO' in both 
cases, and for neutrons with Ej = 41 meV a pyrolytic 
graphite filter was used. For Ej = 115 meV, the experi- 
ment was carried out without a filter in order to maximize 
the neutron flux. Higher order contamination from neu- 
trons with energies above ~ 400 meV is not a concern as 
it results from the high-energy-tail of the thermal neu- 
tron spectrum peaked at <~ 30 meV. The solid lines in 
Figs. 3 and 4 are the result of fits to the 2D Lorentzian 
form 



S(q2D) = 



S(0) 



2 ' 



(5) 



where q^D is the 2D deviation in wave vector from the 
positions of the (1, 0, 1) and (0, 1, 1) rods, convolved with 
the instrumental resolution function of the spectrometer. 

The results so-obtained for the inverse correlation 
length are shown in Fig. 5. These data are consis- 
tent within the errors with the earlier results of Keimer 
et al. [9], but are much more precise and cover a wider 
range of temperatures. The solid line is the predicted 
behavior for the QNLoM in the rcnormalized classical 
regime [6, 7]; this will be discussed below. The results 
for the Lorentzian amplitude 5(0)/£ 2 are shown in Fig. 
6. The four sets of data are normalized to unity over the 
temperature interval 450 K < T < 550 K. 

We now compare the results in Fig. 5 for the correla- 
tion length in La2Cu04 with the predictions of various 
theories. We begin with the results of QMC calculations 
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FIG. 3. Representative energy-integrating two-axis scans 
in La2CuC>4 with Ei = 41 meV and collimations 20' - 10' - S 
- 10'. The solid lines are the result of fits to a 2D Lorentzian 
scattering function Eq. (5) convolved with the resolution 
function of the spectrometer. 

for Eq. (1) with S = 1/2. Because of both advances 
in computational techniques and the implementation of 
finite-size scaling methods, QMC data now exist for £/a 
for the S = 1/2 nn 2DSLQHA for length scales vary- 
ing from 1 to 350,000 lattice constants. QMC results 
of Beard et al. [13] and Kim and Troyer [14] are plot- 
ted in Fig. 7 together with our experimental results in 
La 2 Cu04. The data are plotted in the reduced form £/a 
vs. J/T. It is evident that the QMC and La 2 Cu04 re- 
sults agree in absolute units over the complete tempera- 
ture range (337 K < T < 824 K) or equivalently, length 

scale range (3 ~ £/a ~ 115). Thus, over this range the 
2D spin correlations in La2Cu04 are entirely determined 
by the leading near-neighbor Heisenberg couplings and 
the anisotropic in-plane plus interplanar terms in Eq. (2) 
have no measurable effect to within the uncertainty of our 
experiments. Specifically, the tetragonal-orthorhombic 
structural phase transition at 530 K does not manifest 
itself in the temperature dependence of the correlation 
length. 

We now consider the predictions of various analytical 
theories. A low-temperature theory for the 2DSLQHA 



FIG. 4. Representative energy-integrating two-axis scans 
in La2Cu04 with Ei = 115 meV and collimations 20' - 10' - S 
- 10'. The solid lines are the results of fits to a 2D Lorentzian 
scattering function convolved with the resolution function of 
the spectrometer. 

was formulated by Chakravarty, Halpcrin, and Nelson [6], 
in which they obtained the static and dynamic proper- 
ties of the 2DSLQHA by mapping it onto the 2D quan- 
tum non-linear a model. The 2D QNLcrM is the simplest 
continuum model which reproduces the correct spin- wave 
spectrum and spin-wave interactions of the 2DSLQHA at 
long wavelengths and low energies. First, CHN argued 
that for S > 1/2 the nn 2DSLQHA corresponds to the 
region of the 2D QNLcrM in which the ground state is 
ordered - the renormalized classical regime. Then, CHN 
used perturbativc rcnormalization group arguments to 
derive an expression for the correlation length to two- 
loop order, showing a leading exponential divergence of 
£ versus inverse temperature. Later, Hasenfratz and Nie- 
dermayer [7] employed chiral perturbation theory to cal- 
culate the correlation length more precisely to three-loop 
order. In the RC scaling regime, the correlation length 
is given by 



i 

a 



e c/a 
8 2irp s 



T 



I* 

2 V 2TT Pi 



O 



( T 



\2-rrp. 



(6) 



4 



0.12 



0.10 



0.08 



-0.06 



0.04 



0.02 



Inverse Magnetic Correlation Length 



E, = 14.7 meV 



= 41 meV 
= 90 meV 
= 115 meV 



Renormalized Classical 



I I 



v 



300 400 500 600 700 
Temperature (K) 



800 



FIG. 5. Inverse magnetic correlation length of La2Cu04. 
The solid line is Eq. (6) with J = 135 meV. The Neel and 
structural transition temperatures are indicated by arrows. 

which we refer to as the CHN-HN formula. The parame- 
ters p s and c are the macroscopic T = spin stiffness and 
spin-wave velocity, respectively. For the nearest-neighbor 
2DSLQHA, they are related to the microscopic param- 
eters J, S and the lattice constant a according to p s = 
Z p (S)S 2 J and c = Z c {S)2s/2aSJ . The coefficients Z p (S) 
and Z C (S) are quantum renormalization factors, which 
can be calculated using spin-wave theory (S > 1/2) [22, 
23], series expansion (S = 1/2, 1) [23], and Monte Carlo 
techniques (S = 1/2) [12-14]. For 5 = 1/2, the spin- wave 
approximation gives Z p (l/2) ~ 0.699 and Z c (l/2) ~ 1.18 
[22, 23]. The most precise values for S = 1/2 currently 
available come from the QMC study of Beard et al. [13] 
who find c = 1.657(2)Ja, and p s = 0.1800(5)J and, for 
the T=0 sublattice magnetization, M s = 0.30797(3)/a 2 . 
These correspond to Z c {\/2) = 1.172 and Z Pe {l/2) = 
0.72. The CHN-HN prediction for the Lorentzian ampli- 
tude S(0)/f 2 is 
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It is of interest to compare the QNLctM predic- 
tions with the corresponding predictions of the classical 
spin model for the nn 2DSLHA at low temperatures 
[2,6,24,25] 
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FIG. 6. Lorentzian amplitude, <S(0)/£ 2 versus tempera- 
ture. The data for the different incident neutron energies are 
normalized to unity in the temperature range 450 K ~ T ~ 
550 K. 
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FIG. 7. The logarithm of the reduced magnetic correla- 
tion length £/a versus J/T. The closed circles are data for 
La2Cu04 plotted with J = 135 meV, the open circles are 
data for Sr 2 Cu0 2 Cl2 plotted with J = 125 meV [10], and 
the open squares are the results of the Monte Carlo computer 
simulations [12-14]. The solid line is the theoretical predic- 
tion without adjustable parameters of the 2DQNLcrM for the 
renormalized classical regime Eq. (6). 
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where for classical unit vector spins, p c i = J. For large 
S quantum spins, one finds that the classical limit is ap- 
proached smoothly as a function of S provided that tem- 
perature is measured in units of JS(S+ 1), implying one 
that should take p c i = JS(S + 1) [2]. This choice for 
S = 1/2, gives pd — 0.75J compared with p s = 0.18J. 
The arguments in the exponentials in Eq. (6) and (8) then 
differ by more than a factor of 4 - a very dramatic differ- 
ence between renormalized classical and classical scaling 
behavior. 

An alternative theoretical analysis of the 2DSLQHA 
has been carried out by Cuccoli et al. [26] in which they 
treat quantum fluctuations in a self-consistent Gaussian 
approximation separately from the classical contribution. 
In their approach, which they label the purely-quantum 
self-consistent harmonic approximation (PQSCHA), the 
quantum spin Hamiltonian is rewritten as an effective 
classical Hamiltonian, where the temperature scale is 
renormalized due to quantum fluctuations, and the new 
classical spin length appears as S + 1/2. Defining the 
reduced temperature as t= T/{J(S + 1/2) 2 }, the cor- 
relation length for the 2DSLQHA is then simply given 

by 



£(t)=fci(tcl) with *c 
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Here, £ c ; is the correlation length for the corresponding 
classical 2D square-lattice nn Heisenberg model, which 
is given by Eq. (8) at low temperatures and may be 
obtained from classical spin Monte Carlo calculations 
at higher temperatures [27], and 6* 4 (t) is a temperature 
renormalization parameter. The PQSCHA is most ac- 
curate in the limit where the quantum fluctuations are 
weak, and correspondingly 6* 4 is near unity. This is 
the case for large spin, for example, S=5/2, and in- 
deed this has recently been studied experimentally by 
Leheny et al. [28] who have measured S((]2d) in the 
S = 5/2 2DSLQHA material Rb 2 MnF 4 . They follow a 
field-temperature trajectory which approaches the bicrit- 
ical point in the phase diagram and which accordingly 
should show pure 2D Heisenberg behavior. In this case 
the PQSCHA predicts the correlation length precisely 
with no adjustable parameters over the inverse temper- 
ature range 0.5 < p c i/T < 2 or equivalently, the length 



scale range 1< £/a ~ 100. We note from Eq. (8) and (9) 
that for the classical model T always appears scaled by 
Pd. Thus the quantum effects in the PQSCHA can be 
thought of simply as a temperature dependent renormal- 
ization of pd, that is, pd — > 9 4 (t)p c i. 

Finally, for the QNLctM there may be a crossover from 
renormalized classical to quantum critical behavior with 
increasing temperature [6]. In the QC regime heuristi- 
cally one expects 

c/a 



£/a = 0.8- 
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T-T QC 

with Tqc > adjustable [2,6,10]. We emphasize that 
this anticipated crossover is a property of the QNLctM 
and it may or may not occur for quantum spins on a 2D 
lattice. 

The solid line in Fig. 5 is the QNLctM-RC prediction, 
Eq. (6) with c and p s from Beard et al. [13] As ob- 
served previously for Sr 2 Cu02Cl2 [10], Eq. (6) describes 
the measured correlation length of La2Cu04 extremely 
well without adjustable parameters over the temperature 
range 337 K < T < 824 K, or equivalently, the length 
scale range, ~ 3 ~ £/a^ 110. All of the data for £/a 
from each of quantum Monte Carlo, Sr 2 Cu02Cl2 and 
La2Cu04 together with Eq. (6) are plotted in the uni- 
versal form £/a vs. J/T in Fig. 7. The evident universal 
behavior is, of course, both pleasing and reassuring. The 
good agreement of all of the experimental and numerical 
results with the low temperature QNLctM- RC predictions 
down to very small values of J/T at first appears to be 
quite puzzling. The QMC study, Ref. [13], suggests that 
this agreement is, at least in part, accidental. Specifi- 
cally, in crossing over from the low temperature contin- 
uum QNLctM to the discrete lattice S=l/2 Heisenberg 
model the higher order terms in Eq. (6) conspire such 
that over the measured temperature range the deviation 
of £/a from Eq. (6) is never more than 15% which is well 
within the experimental error. 

We now focus on the high temperature behavior in 
La2Cu04. We show in Fig. 8 the La2Cu04 correlation 
length data together with the predictions from QNLctM- 
RC (Eq. (6)), QNLctM-QC Eq. (11), high temperature 
series expansion [29] and the PQSCHA which involves 
Eq. (10) combined with results of classical Monte Carlo 
simulations [27]. As observed previously for Sr2Cu02Cl2 
[10] as well as for both S = 1/2 2DSLQHA QMC cal- 
culations [14] and high temperature series expansion re- 
sults [29], the QNLctM-QC prediction, Eq. (11), disagrees 
strongly with the experimental results in La2CuC>4. This 
is, perhaps, not surprising given the extremely short 
length scales at the relevant temperatures. Specifi- 
cally, at these short distances, the continuum QNLctM 
approach which underlies the possible QC behavior is 
probably no longer valid. By contrast the PQSCHA 
which corresponds to classical scaling for the pure 2D 
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FIG. 8. The logarithm of the reduced magnetic correla- 
tion length £/a versus J/T compared with the predictions of 
various theories including renormalized classical behavior, Eq. 

(6) , quantum critical behavior, Eq. (11), the PQSCHA, Eq. 
(10), and high temperature series expansion [29]. 

Heisenberg model [2,30] agrees reasonably well in abso- 
lute units with no adjustable parameters for length scales 
up to about £/a <~ 15. As expected, the PQSCHA breaks 
down at lower temperatures and larger length scales. 
Thus, if there is a crossover in the correlation length, it 
is from classical scaling to renormalized classical scaling 
with decreasing temperature. Clearly, it is very impor- 
tant that theory for this crossover from the high temper- 
ature PQSCHA classical scaling regime to the low tem- 
perature QNLcrM-RC regime be developed. 

Finally, we discuss the behavior of the structure factor 
S(0). The leading divergence of 5(0) is determined by 
£ 2 . This is confirmed by the results for 5(0) in La2Cu04 
displayed in Fig. 6 which shows that 5(0)/£ 2 is approx- 
imately constant over the complete temperature range. 
This is equivalent to the statement that for the 2D S=l/2 
QHA the critical exponent 772 = 0. We note that in 
3D r/ 3 ~ 0.04, [31] that is, 5(0) - £ 1 - 96 whereas for 
the ID S=l/2 QHA one has the remarkable result that 
5(0) - {Inif/ 2 [32], which implies m = 2. 

The temperature dependent correction factors (c.f. Eq. 

(7) ) beyond the leading £ 2 divergence are problematic. 
Specifically, Greven et al. [10] find in their measurements 
in Sr 2 Cu02Cl2 that over the length scale 5~ £/a~200, 
5(0)/£ 2 is independent of temperature to within the er- 



rors. By contrast, QMC [12-14] and high temperature se- 
ries expansion [29,30] studies of the S=l/2 nn 2DSLQHA 
find 5(0)/£ 2 ~ T 2 over about the same range of length 
scales. In the S = 5/2 2DSLQHA material Rb 2 MnF 4 , 
Leheny et al. [28] find a clear crossover at £/a~4 from 
5(0)/£ 2 ^T 2 behavior to a much weaker dependence of 
5(0) fi 2 on T. The data for 5(0) /£ 2 in La 2 Cu0 4 shown in 
Fig. 5 are clearly inconsistent with T 2 behavior over the 
complete temperature range but would allow a gradual 
crossover as found in Rb 2 MnF4. This lack of universality 
in 5 (0) / £ 2 seems surprising given the robust universality 
of the behavior for £/a (Fig. 7). Of course, departures 
from the low temperature QNLcrM-RC behavior may oc- 
cur at different temperatures for different quantities. It 
is also possible that the terms in Eq. (2) beyond the 
nearest neighbor Heisenberg coupling will effect 5(0) /£ 2 
more than they effect £ itself. 

In summary, we have carried out a neutron scat- 
tering study of the instantaneous spin-spin correlations 
in La 2 Cu04 (Tjy=325 K) over the temperature range 
337 K to 824 K. Incident neutron energies varying from 
14.7 meV to 115 meV have been employed in order to 
guarantee that the energy integration is carried out prop- 
erly. The results so-obtained for the spin correlation 
length as a function of temperature when expressed in 
reduced units agree quantitatively both with previous 
results for the 2D tetragonal material Sr 2 Cu0 2 Cl 2 and 
with quantum Monte Carlo results for the nearest neigh- 
bor square lattice S = 1/2 Heisenberg model. All of the 
experimental results for the correlation length are well 
described without any adjustable parameters by the be- 
havior predicted for the quantum non-linear sigma model 
in the low temperature renormalized classical regime. 
The structure factor, on the other hand, deviates subtly 
from the predicted low temperature behavior although 
the leading £ 2 behavior is confirmed. The correlation 
length data at high temperature agree reasonably well 
with predictions of the PQSCHA which corresponds to 
classical scaling with quantum corrections for the 2D 
Heisenberg model. We therefore hypothesize that in 
La 2 Cu04 there is a gradual crossover from renormalized 
classical to classical scaling with increasing temperature. 
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